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In this paper we propose a model to reproduce superconductive and normal properties of the iron 
pnictide LiFeAs in the framework of the four-band s± wave Eliashberg theory. A confirmation of the 
multiband nature of the system rises from the experimental measurements of the superconductive 
gaps and resistivity as function of temperature. We found that the most plausible mechanism 
is the antiferromagnetic spin fluctuation and the estimated values of the total antiferromagnetic 
spin fluctuation coupling constant in the superconductive and normal state are Atot = 2.00 and 
Atot.tr = 0.77. 

PACS numbers: 74.70.Xa, 74.25. F, 74.20. Mn, 74.20.-z 



Recent ARPES measurements of iron superconductor 
LiFeAs report four slightly anisotropic gaps Their 
isotropic values at 8 K are given by Ai = 5.0 meV, A2 = 
2.6 meV, A3 = 3.6 meV, A4 = 2.9 meV and the critical 
temperature for this compound is T c — 18 K [2|. 

In an other work Q we disregarded the anisotropic 
part of the gap values and we tried to reproduce the ex- 
perimental data in the framework of s± wave multiband 
Eliashberg theory. At first, we calculated [3-8] the Fermi 
surface, depicted in FIG. [TJ Five different sheets are 
present, with two electron pockets centered near the M- 
point of the Brillouin zone and three hole pockets around 
the r-point. The 5-th sheet can be disregarded thanks 
to its low density of states and size Q as can be seen 
in TABLE U In this way a four-band s-wave Eliashberg 
model [9|, llfj can be used and eight coupled equations 
for the gaps Aj(m? n ) and the renormalization functions 
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TABLE I: Fermi Surface resolved Kohn Sham properties: the 
Fermi density of states N(0) is given in states/spin/eV and 
plasma frequencies u p in eV. ab label the in-plane and c for the 
out-of-plane direction of the diagonals of the plasma tensor 



Ziiiuin) have to be solved. If i is the band index (that 
ranges between 1 and 4) and uj n are the Matsubara fre- 
quencies, the imaginary-axis equations are: 

+ £[ r ij + Af]^(^«), (i) 
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FIG. 1: Fermi surface of LiFeAs Q. 



where Tjj and Tff are the non magnetic and magnetic 
impurity scattering rates, Q(u> c — \uj m \) is the Heaviside 
function and w c is a cutoff energy. Moreover, h*j(ui c ) 
are the elements of the 4x4 Coulomb pseudopoten- 
tial matrix and N^(iuj m ) = Aj(iuj m )/ ^Jlj^ + A|(iw m ), 

Nf(iu m ) = u m / Ju^ + A? {iw m ). Finally, 
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A^(iw„, iuj m ) = A?j(iuj n , iuj m ) - A*f(iui n , iu> m ). 

Here the superscripts sf and ph mean "antiferromagnetic 
spin fluctuations" and "phonons" , respectively. In par- 
ticular, 

ij =*J o dnn K - a , m) 2 + n 2 . 

and the electron-boson coupling constants are defined as 
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TABLE II: In the first row are shown the experimental data. 
The second row concerns the pure intraband case (Am = 0.0) 
while the third one concerns the case with a large intraband 
term (An = 2.1). The critical temperatures are given in K 
and the gap values in meV. 
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The solution of eqs.([T|) and @ requires a huge number 
of input parameters, then drastic approximations, com- 
patible with the goal of reproducing the essential physics 
of the problem, are necessary to make the model solvable. 
As for many other pnictides we assumed that [3j: i) the 
total electron-phonon coupling constant is small [111: ii ) 
spin fluctuations mainly provide interband coupling [l2j . 
This means that we can set A^ — A?- — 0, p*i{oJ c ) — 
/i*j(uj c ) — 0, i.e. the electron-phonon coupling constant 
and the Coulomb pseudopotential in first approximation 
compensate each other and X S J — (only interband SF 
coupling) However, within these assumptions, we 

were not able to reproduce the observed gap values, and 
in particular the high value of Ai. In order to solve this 
problem it is necessary to introduce an intraband cou- 
pling in the first band (An ^ 0). 

The final matrix of the electron-boson coupling constants 
becomes 
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\ A41 = Ai4^i4 A 42 = A24^24 
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where i/y = Ni(0)/Nj(0) and Ni(0) is the normal den- 
sity of states at the Fermi level for the i-th band (i = 
l,2,3,-£h_ We choose spectral functions with Lorentzian 
shape 



a%Fij (fi) = C l3 {L{9, + Qij , Y iS ) - L(fi 



tti^Yij)} (5) 



5- and Cij are 



where L{Q ± fly.Yy) = (0±0<j)a+y< j 
normalization constants, necessary to obtain the proper 
values of A^ while and YL- are the peak energies 
and half-widths of the Lorentzian functions, respectively 
121. In all the calculations we set fL, = f2„ s „ = fin 



meV [3, and Y l} = Y t s / = n s J /2 (ij. The cut-off 
energy is ui c = 18 Qq and the maximum quasiparticle 



energy is w„ 
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Bandstructure calculations 



(see TABLE fl} provide information about the factors 
Vij that enter the definition of A,j. In the end the 
model contains five free parameters: The coupling 



constants A13, A23, A14, A24 and An. First of all we 
solved the imaginary-axis Eliashberg equations ([1]) and 
@ (actually we continued them analytically on the 
real-axis by using the Pade approximant technique) and 
we fixed the free parameters in order to reproduce the 
gap values at low temperature. The large number of 
free parameters (five) may suggest that it is possible to 
find different sets that produce the same results. On 
the contrary, as a matter of fact, the predominantly 
interband character of the model drastically reduce 
the number of possible choices. At this point there 
are no more free parameters. We can calculate the 
critical temperature and it turns out to be very close 
to experimental one @: T c c alc = 18.6 K. In TABLE |TT] 
the obtained results are summarized. The problem of 
this model is the necessity of a so large intraband term 
An in order to give a physical interpretation of the 
experimental data 0]. 

Regarding the normal state [l5| , the resistivity satura- 
tion at high temperature [l6[ suggests that the presence 
of several sheets in the Fermi surface also affects the nor- 
mal state transport properties. 

First of all we noticed (see FIG. [2]) that at low tem- 
perature p{T) oc T 2 and this could indicate that a non- 
phononic mechanism plays a relevant role in the physics 
of this compound [17 1 . 

To begin with, we tried to fit the data within a one-band 
model 3, 19[ (see eq. ([6]) with i = 1) where the phonon 
spectra has been taken from ref. [20( and the plasma 
energy has been obtained by first principle calculation 
(see TABLE fl}. The transport coupling constant and the 
value of the impurities are considered as free parame- 
ters. The obtained values are reported in TABLE IlIIl in 
particular At r ,tot = 0.32 which is in agreement with the 
calculated value of the trasport electron-phonon coupling 
constant 2l|. However, as can be seen in FIG.[5J within 
a one-band model (black dashed line) the experimental 
data cannot be reproduced. 

Phenomenological model [22} proposed to explain satura- 
tion at high temperature generally assume the presence 
of parallel conductivity channels where one of them has a 
strong temperature dependence and another one is char- 
acterized by a temperature-independent contribution. 
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TABLE III: The first and second rows concern the phonon 
case while the third one concerns the case of the antiferro- 
magnetic spin fluctuation spectral function. The 71 and Qo 
are given in meV. 
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FIG. 2: Temperature dependence of resistivity in LiFeAs. Ex- 
perimental data and calculated fits are reported. The black 
dashed line comes from a single-band model. Within a four- 
band model two different cases have been considered, one is 
obtained with the phononic spectra (green dash-dotted line) 
and one with the antiferromagnetic spin fluctuation spectra 
(red solid line). The inset shows the two normalized spectral 
function that have been used, the phonon spectra (black solid 
line) and the antiferromagnetic spin fluctuation spectra (red 
dashed line) 



In the wake of the model proposed for the superconduct- 
ing state, we propose a multiband model 23j, |24j for ana- 
lyzing the resistivity data. We will examine two possible 
mechanism responsible of resistivity: Phonons and anti- 
ferromagnetic spin fluctuations. The theoretical expres 



sion of the resistivity as function of temperature 
is given by the equation: 
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where Wp^j is the bare plasma frequency of the i-band 
and 
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here 71 = X^=i + Pjjf is the sum of the inter- and 
intra-band non magnetic and magnetic impurity scatter- 
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FIG. 3: Magnification of the previous figure. Resistivity at 
low temperature. 



ing rates present in the Eliashberg equations and 



x> 2 

.7=1 



tr,ij 



(8) 



where a^ r (f2).Ftx,ij (Q) are the transport spectral functions 
related to the Eliashberg functions [3| • 
If a normalized transport spectral function a 2 T (Q)F{ T ;(f2) 
is defined, then c4(fi)F tr)ij (fi) = Atr,«a? r (n)F^Jj(n) 
where the coupling constants are defined as for the stan- 
dard Eliashberg functions. 

In order to build a model as simple as possible, we 
chose all the normalized transport spectral functions to 
be equal, then a\ Y (Sl)F' tx ;(fi) = At r ,ia? r (fi)F t ' r (fi) where 

It has been shown that, at least for iron pnictides, this 
model can have a theoretical support [24| depending on 
the electronic structure of the compound. The basic idea, 
based on ARPES and de Haas-van- Alphen data, is that 
the transport is drown mainly by the electronic bands 
and that the hole bands have a weaker mobility [2^ . 
Then the impurities are mostly present in the hole bands 
and 712 3> 73,4, while the transport coupling is much 
higher in bands 3 e 4 and this means that, at least as a 
first approximation, Ai and A2 can be fixed to be zero. 
In this way we will have two contributions almost tem- 
perature independent and two which change the slope of 
the resistivity with the temperature [24 1 . 
Let us start with the phononic case. For simplicity we 
considered all the spectral function to be proportional 
to the phonon spectra used also in the previous fit [20j ]. 
As mentioned above the transport spectral functions are 
similar to the standard Eliashberg functions. The main 
difference is the behavior for Q, — ► [3] , where the trans- 
port function behaves like Jl 4 instead of tt 2 as in the su- 
perconducting state. So the condition a 2 r (ri)F tr (ri) oc il 4 
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has been imposed in the range < ft < fceTb/lO and 
then 

a4(Q)F t ' r (Q) = 6iO 4 i?(fc B T D /10 - fi) 

+ Ci al.(n)Fl' I (n)$(n - fc B T D /io), 

where = 240 K is the Debye temperature [26|, the 
constant 6; and c; are fixed by imposing the continuity in 
fcB?b/10, and the normalization to 1 and a 2 
is proportional to electron-phonon spectral function 
al(ft)F{ r (n) is shown in the inset of FIG. H 
All the plasma frequencies are fixed by first principle cal- 
culations (see TABLE I) and the coupling constants are 
considered as free parameters as well as the impurities pa- 
rameters. The best fit is obtained with At r ,tot = 0.14, as 
reported in TABLE [TTT1 which is in agreement with the 
hypothesis that the phonon coupling in LiFeAs is very 
weak and the value of XtrA almost does not influence the 
final result. However the experimental data are not per- 
fectly reproduced, as can be seen by looking the green 
dash-dotted curve in FIG. [3] and in FIG. |3l moreover a 
huge quantity of impurity has been necessary to obtain 
this theoretical curve and this is not in agreement with 



of the spectral function and the total transport coupling 
constant have been extracted. 



16] 



the good quality of single crystal 
Then we considered the case of antiferromagnetic spin 
fluctuations. Now for O — ^ the transport function be- 
haves like J7 3 instead of f2 as in the superconducting state. 
So the condition a^ I (^l)F tr (^l) oc O 3 has been imposed 
in the range < Q < fl /10 and then a$ T (Q)Fl T (Q) = 
bifl 3 ^(n /l0 Cial(n)F^(Q)^(Q - n /10) and the 

constant b\ and c; are fixed in the same way as before. 
We choose as a| r (Q)F t "(Q) the theoretical antiferromag- 
netic spin fluctuation function in the normal state 27 j 



alF"(n) oc 



n n 



rtf(fi-n ), 



(9) 



n 2 + n 2 

where fio is a free parameter: from the fit of experimen- 
tal data we obtain fio = 47 meV. 

Also in this case the value of the free parameters are re- 
ported in TABLE [TTT] and FIG. H depicts the obtained 
results with the red solid line as well as the spectral 
function (in the inset). The curve obtained by using the 
spin fluctuation spectra better reproduce the experimen- 
tal with a total coupling given by A tr tot = 0.77 consis- 
tent with expectations, indeed is smaller than the value 
in the superconducting state. Moreover the parameters 
seems to better represent the LiFeAs sample: this is a 
stoichiometric compound and the data have been taken 
from measurements on a single crystal sample, then the 
presence of a huge amount of impurities is not supported. 
Of course we have done a draft simplification because the 
more plausible situation is the coexistence of two mech- 
anisms but certainly the antiferromagnetic spin fluctu- 
ations constitute the main mechanism. In conclusion 
we can say that in this compound the antiferromagnetic 
spin fluctuations play an important role also in the nor- 
mal state, moreover information about the energy peak 
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